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O Abstract 
(N 

O First measurements of azimuthal asymmetries in hadron-pair production in deep-inelastic scattering 

1^ of muons on transversely polarised ^LiD (deuteron) and NH3 (proton) targets are presented. The data 

Mh were taken in the years 2002-2004 and 2007 with the COMPASS spectrometer using a muon beam of 

00 160 GeV/c at the CERN SPS. The asymmetries provide access to the transversity distribution func- 

tions, without involving the Collins effect as in single hadron production. The sizeable asymmetries 
measured on the NH3 target indicate non-vanishing M-quark transversity and two-hadron interference 
^ fragmentation functions. The small asymmetries measured on the ^LiD target can be interpreted as 

D indication for a cancellation of u- and t/-quark trans versities. 
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1 Introduction 

The quark content of the nucleon at twist-two level can be characterised by three independent distri- 
bution functions (DFs) for each quark flavour |[T1. If the quarks are collinear with the parent nucleon, 
i.e. quarks have no intrinsic transverse momentum kr, or after integration over kr, these three distri- 
butions exhaust the information on the partonic structure of the nucleon lHJ^. Two of these functions 
have been measured with good accuracy. First, the unpolarised distribution ff{x) describes the number 
density of quarks with flavour q carrying a fraction x of the longitudinal momentum of the fast moving 
parent nucleon. Secondly, the helicity distribution g\{x) describes the difference between the number 
densities of quarks with helicity parallel and antiparallel to the spin of the nucleon that is longitudinally 
polarised with respect to the nucleon momentum. However, up to ten years ago nothing was known 
about the third function, the transverse spin distribution often referred to as transversity, which 

describes the difference between the number densities of quarks with their spins oriented parallel and 
antiparallel to the spin of a transversely polarised nucleon. This distribution is difficult to measure, since 
it is related to soft processes correlating quarks with opposite chirality, making it a chiral-odd function 
m . As a result, transversity can only be accessed through observables in which it appears coupled with 
a second chiral-odd object in order to conserve chirality, i.e. not in inclusive deep inelastic scattering 
(DIS) at leading twist. For a recent review see Ref. |6|. 

The only results on transversity so far come from semi-inclusive deep inelastic scattering (SIDIS) reac- 
tions. Here, the chiral-odd partners of the transversity distribution function are fragmentation functions 
(FFs), which describe the spin-dependent hadronisation of a transversely polarised quark q into hadrons. 
Already in 1993 it was suggested by Collins [7 | that transversity could be accessed in SIDIS due to a 
hadronisation mechanism, in which the transverse spin of the struck quark is correlated to the transverse 
momentum pr of the produced unpolarised hadron with respect to the quark momentum. The resulting 
azimuthal asymmetry in the distribution of the final state hadrons is then proportional to the convolution 
of this Collins FF Hy{z,pT) and the transversity DF h\{x,kT). Here, convolution integrals of quark and 
hadron transverse momenta have to be taken into account [3lil4l[8|. 

Such asymmetries were measured by HERMES using a transversely polarised pure hydrogen target ||9l 
and by COMRA^SS using transversely polarised ^LiD (deuteron) liTOl and NH3 (proton) targets ifTTI . No 
significant signals were found on the ^LiD target. Sizeable Collins asymmetries were observed by HER- 
MES, and by COMR\SS using the polarised NH3 target, which implies both transversity DF and Collins 
FF to be different from zero. The latter was also independently measured at the KEK e^e^ collider by 
the BELLE experiment and established to be sizeable [[1211131 . All these results were well described by 
a global fit IIT4IIT5I . which allowed for a first extraction of the u- and <i-quark transversity DFs. These 
results contain convolutions of transverse-momentum-dependent functions, where the assumed shape of 
the transverse momentum distributions leads to a model dependence of the extracted results. Moreover, 
BELLE results as compared to COMR^SS and HERMES results were obtained at a very different scale, 
« 100(GeV/c)2 compared to {Q^) w 3.2(GeV/c)2 and {Q^) w 2.4(GeV/c)2, respectively. The nec- 
cessary evolution is quite complex for transverse-momentum-dependent functions |l8lfT6l, requiring 
the use of CoUins-Soper factorisation lITTl [181 and thus again the knowledge of the transverse momen- 
tum distributions. In the global analysis of Refs. |[T4l[T5l . a gaussian ansatz for the transverse momentum 
distributions was used but only collinear evolution was taken into account, leading to a possible overes- 
timation of h\ [8|. 

As an alternative to study the Collins fragmentation mechanism in one-hadron SIDIS it was proposed 
to measure two-hadron SIDIS, Ip^ I'h^h^X with both hadrons produced in the current fragmentation 
region ||T9l - [22ll . In this reaction appears a new chiral odd fragmentation function, the Interference Frag- 
mentation Function (IFF) H^, which describes the fragmentation of a transversely polarised quark into 
a pair of unpolarised hadrons. The transverse polarisation of the fragmenting quark is correlated with 
the relative momentum of the two hadrons, which gives rise to an azimuthal asymmetry with respect to 
the virtual-photon direction and the lepton scattering plane. The additional degrees of freedom allow for 
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2 THEORETICAL FRAMEWORK 



an integration over the transverse momenta of the final state hadrons, leaving only the relative momen- 
tum of the two hadrons. This avoids the complexity of transverse-momentum-dependent convolution 
integrals and the analysis can be performed using collinear factorisation fSl [T3 . Here, the evolution 
equations are known at next-to-leading order |[23l . so that results from e^e^ scattering and SIDIS can 
be correctly connected, making it the theoretically cleanest way to extract transversity using presently 
existing facilities [8J. The properties of interference fragmentation functions are described in detail in 
Refs. UMllEHlIl. 

First evidence for azimuthal asymmetries in leptoproduction of tt+tt^ pairs was published by HERMES, 
using a transversely polarised hydrogen target |28|. The interference fragmentation functions were mea- 
sured in e^e^ reactions by BELLE [29 1. These measurements indicate a sizeable M-quark transversity 
distribution and non- vanishing interference fragmentation functions ll30ll . In this Letter, the first mea- 
surement of two-hadron azimuthal asymmetries using a transversely polarised ^LiD (deuteron) target is 
presented as well as results from a NH3 (proton) target. Due to the large acceptance of the COMRA^SS 
spectrometer and the large beam momentum of 160GeV/c, results with high statistics were obtained 
covering a large kinematic range in x and Mhh, the invariant mass of the hadron pair. 



2 Theoretical Framework 

At leading twist and after integration over the hadron transverse momenta, the cross section of semi- 
inclusive two-hadron leptoproduction on a transversely polarised target is given as a sum of a spin inde- 
pendent and a spin dependent part |[25l 



J cos Q dM^i^ d(^R dz dx dy d(j)s ^nQ^y \ 2 

xY^elf,{x)Di^,{z,MLcosd) 



1-3'+^^ (1) 



d^o^J■ _ a} 



S^{\-y) (2) 



d cos Q dM\y^ d <pR dz dx dy d^s ^.TiQ^y 

Here, the sums run over all quark flavours q, and pi and p2 denote the momenta of the two hadrons of 
the pair. The first index (U) indicates an unpolarised beam and the second, (U) or (T), an unpolarised 
and transversely polarised target, respectively. Note that the contribution from a longitudinally polarised 
beam and a transversely polarised target, a^j-, is suppressed by 1/Q ||3TI . The fine-structure constant is 
denoted by a, y is the fraction of the muon's energy lost in the nucleon rest frame, Di^(z,M^^,cos0) is 
the unpolarised two-hadron fragmentation function for a quark of flavour q and zi, Zi are the fractions 
of the virtual-photon energy carried by these two hadrons, with z = Zi+Zi- The symbol denotes the 
component of the target spin 5 perpendicular to the virtual-photon direction, and Q is the polar angle 
of one of the hadrons - commonly the positive one - in the two-hadron rest frame with respect to the 
two-hadron boost axis. The azimuthal angle (^rs is defined according to Fig. [T] as 



= ^R-^s' = h + ^s-^, (3) 

where is the azimuthal angle of the initial nucleon spin and 0y is the azimuthal angle of the spin of 
the fragmenting quark, with (ps' = 7i — (ps- The azimuthal angle (pR is defined by 

^ - (^"^^-^ arccosf^i^^V (4) 



\{qxl)-R\ \ \qxl\\qxR 
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Fig. 1 : Definition of the azimuthal angles 0r and 05 for two-hadron production in deep inelastic scatter- 
ing, where /, q and p,- are the 3-momenta of beam, scattered muon, virtual photon and hadrons. Note 
that the azimuthal plane is defined by the directions of the relative hadron momentum and the virtual 
photon. 



where I is the incoming lepton momentum, q the virtual-photon momentum and R the relative hadron 
momentum ll20l [32l given by 

R = , =: S2P1 - SiP2 ■ (5) 

Zl +Z2 

The number Nh+h- of pairs of oppositely charged hadrons produced on a transversely polarised target 
can be written as 

Nh+h- {x,y,z,Mhh , COS , (j)Rs) cc Ouu (^±f{x,y)PT Dnn A^r*^""' sin sin ) , (6) 



omitting luminosity and detector acceptance. Here, Pt is the magnitude of the transverse target polari- 
sation and Dnn{y) = ^_y^y2/7 the transverse-spin-transfer coefficient, while f{x,y) is the target dilution 
factor that depends on kinematics. It is given by the number-weighted ratio of the total cross section for 
scattering on protons or deuterons to that for scattering on all nuclei in the target. It increases at large x 
due to the reduced cross section for heavy targets. The dependence of the dilution factor on the hadron 
transverse momenta appears to be weak in the v range of the COMPASS experiment. Dilution due to 
radiative events on protons or deuterons is taken into account by the ratio of the one-photon exchange 
cross section to the total cross section. For ^LiD, / includes a correction for the relative polarisation 
of deuterons bound in ^Li with respect to free deuterons. For ^'*NH3, / contains a correction for the 
polarisation of the small ^^N admixture. 



The asymmetry 



" 2Mhh 1,^2 -/fCx) •Di,,(z,M2„ cos 0) 

is proportional to the product of the transversity distribution function h\ {x) and the polarised two-hadron 
interference fragmentation function H^^{z,Mf^jj, cos 6), summed over the quark flavours q with charge 
eg. 

It is convenient 1251 to expand both the polarised and unpolarised two-hadron fragmentation functions 
in terms of Legendre polynomials in cos0. For the invariant mass range typically covered by SIDIS 



\pi-p2\L,el-h1{x)-H<^{z,Ml,cosd) 
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3 EXPERIMENTAL DATA 



experiments, M/,/, < 1.5 GeV/c^, to a good approximation only relative s and p partial waves of the two- 
hadron system contribute to the cross section, yielding 1 271 

Di{z,M^^,cose) ~ Di^ooiz,Mlh) + coseDi_oi{z,Mhh) 

+ i(3cos2 0-l)Di,,/(z,M,\) (8) 

and 

H<{zMh,cose)c^H<Jz,ML) + ^oseH<,iz,ML), (9) 

respectively. The term cos dDi oi{z,Mf^f^) describes the interference between an unpolarised hadron 
pair (denoted o) in s-wave and a longitudinally polarised pair (denoted /) in p-wave, while the term 
^toA^^^hh) from the interference between an unpolarised hadron pair in s-wave and a transversely 
polarised pair (denoted t) in p-wave. The term cos 6 H^i^{z,Mf^f^) indicates interference between longi- 
tudinally and transversely polarised pairs in a relative p-wave, while the term (3cos^ 6 — l)Di //(z,M^^) 
indicates interference between longitudinally polarised pairs in a relative p-wave. The term Di^oo{z,M^i^) 
represents an unpolarised state of the hadron pair and can have contributions from s- and p-waves but not 
from the interference between both. 

3 Experimental Data 

The analysis presented in this Letter is performed using data taken in the years 2002-2004 and 2007 
with the COMPASS spectrometer [33] by scattering positive muons of 160GeV/c from the CERN SPS 
off transversely polarised solid state ^LiD and NH3 targets, respectively. The beam muons originating 
from 71+ and decays are naturally polarised with an average longitudinal polarisation of about 0.8 
with a relative uncertainty of 5 %. For ^LiD, the average dilution factor is (/) ~ 0.38 and the average 
polarisation is {Pt) ~ 0.47, while for NH3 it is (/) ~ 0.15 and (Pt) ~ 0.83, respectively. The target con- 
sists of cylindrical cells in a row, which can be independently polarised. In 2002-2004, two cells were 
used, each 60 cm long and 3 cm in diameter. The direction of polarisation in the downstream cell was 
chosen oppositely to the one in the upstream cell. In 2007, the target consisted of three cylindrical cells, 
with 4 cm diameter. The middle cell was 60 cm long and the two outer cells 30 cm each. The direction 
of polarisation in the middle cell was opposite to the one in the outer cells. For the analysis the central 
cell is divided into two parts, providing four data samples with two different orientations of polarisation. 
Both target configurations allow for a simultaneous measurement of azimuthal asymmetries for both tar- 
get spin states to compensate flux dependent systematic uncertainties. Furthermore, the polarisation was 
destroyed and built up in reversed direction every four to five days. This compensates acceptance effects 
created by the dipole field necessary to sustain transverse target polarisation. 

For the analysis events with incoming and outgoing muons and at least two reconstructed hadrons from 
the primary vertex inside the target cells are selected. Equal flux through the whole target is achieved by 
requiring that the extrapolated beam track crosses all cells. In order to select events in the DIS regime, 
cuts are applied on the squared four-momentum transfer, > I (GeV/c)^, and on the invariant mass of 
the final hadronic state, W > 5 GeV/c^. Furthermore, the fractional energy transfer by the virtual photon 
is required to be y > 0. 1 and y < 0.9 to remove events with poorly reconstructed virtual photon energy 
and events with large radiative corrections, respectively. 

The hadron pair sample consists of all combinations of oppositely charged hadrons originating from the 
reaction vertex. Hadrons produced in the current fragmentation region are selected requiring z > 0.1 for 
the fractional energy of each hadron and xp > 0.1. The Feynman variable Xf is defined in the centre- 
of-mass frame of virtual photon and target nucleon as the longitudinal momentum of the hadron with 
respect to the virtual photon divided by the total available longitudinal momentum. The suppression of 
exclusive hadron pairs lITTll is accomplished by requiring the missing mass to be Mx > 2.4GeV/c^. As 
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the azimuthal angle is only defined for non-coUinear vectors R and q, a minimum value is required on 
the component of i? perpendicular to q, \R±\ > 0.07 GeV/c. After all cuts, 5.8 x 10^ h^h^ combinations 
for the ^LiD target and 10.9 x 10^ h^h^ combinations for the NH3 target remain. Figure [2] shows the 
invariant mass distributions of the two-hadron system for both targets, always assuming the pion mass 
for each hadron. A cut of M/,/, < 1.5 GeV/c^ is applied in order to justify the restriction to relative s and 
p waves given in Eqs. ([8]l and (|9]). For further details on the event selection and on the analysis we refer 
to Ref. Il34l . 

In Refs. Il2ni22l|25l . it was proposed to measure two-hadron production integrated over the angle 6. This 
has the advantage that in the resulting expression for the cross sections only the fragmentation function 
H^gf{z,Mf^f^) appears, provided that the experimental acceptance is symmetric in cos 6. 
Figure [3] shows the sin0, cos0 and cos^ 6 distributions for the NH3 target. The corresponding distri- 
butions for the ^LiD target show exactly the same shape. In the COMR\SS acceptance, the opening 
angle 6 peaks close to 7i/2 with (sin0) = 0.94 (Fig. [3j left). In the analysis, we extract the product 
A = {A^^Y^'^^ sin0), integrated over the angle 6. The cosd distribution is symmetric around zero (Fig. 3 
centre) with (cos0) = 0.01. Therefore, the contribution to the asymmetry arising from H^i^{z,Mfl^ 
in Eq. (|9]l is expected to be rather small, so that the result is mainly sensitive to fi^oti^^^lih)- 
cos^ B distribution (Fig. [3] right) shows a mean value of (cos^ d) =0.11. Therefore, the contribution of 
D\ji{z,Mfjfj) in Eq. ([8) does not vanish but contributes with a weight of about 16 % to the unpolarised 
cross section. This leads to a dilution of the asymmetry signal which has to be taken into account when 
extracting the transversity distributions from the data. For this purpose, the mean values of all three 
distributions for individual kinematic bins can be found on HEPDATA ll35l. 



4 Asymmetry Extraction 

The asymmetry A = {Pt)A is evaluated in kinematic bins of x, z or M/,/j, while always integrating over 
the other two variables. As estimator an extended unbinned maximum likelihood function in and 
is used: 




Fig. 2: Invariant mass distributions of the final h^h^ hadron samples for the NH3 target (upper curve) 
and the ^LiD target (lower curve). The cut M/j/, < 1 .5 GeV/c^ is indicated. In both cases the K^, p and f\ 
resonances are clearly visible. 
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4 ASYMMETRY EXTRACTION 




Fig. 3: Distributions in sin 0, cos 6 and cos^ 6 for the h^h sample from the transversely polarised NH3 
target. 



(10) 



where the probabihty density function ^5; af , A) = af ^5) • (1 ± / • D„„ • A • sin (I)rs) is nor- 

malised to the estimated number of hadron pairs if = J J d<^Rd<^sP^{'^R^ '^s\<^f }^)- Here it denotes the 
sign of the target polarisation and af{<^R,<^s) represents the acceptance seen by particles produced in 
target cell /, including the unpolarised cross section and the respective luminosities. The outer product 
corresponds to the target cells, riceu = 2 in case of the ^LiD target and riceu = 4 in case of the NH3 target, 
while the inner products correspond to the two data samples of each cell / acquired with target spin up 
and target spin down, respectively. The contributions are weighted with powers of N/Nf to account for 
unbalanced statistics, where A^^ is the average number of pairs per sample and A^^ is the number of pairs 
with spin up or spin down for each target cell /. The weighting makes the resulting asymmetries less 
sensitive to acceptance and reduces possible false asymmetries. Monte Carlo studies have shown that the 
functional form of the acceptance has a negligible effect on the extracted asymmetries. Hence constant 
values af are used, which account for the different luminosities in the two periods with opposite target 
polarisation. Additionally, each pair is weighted by the corresponding target dilution factor f{x,y) and 
the transverse spin-transfer-coefficient D„„{y). The magnitude of the target polarisation cannot be used 
as a weight in the fit since this could bias the results in case of unbalanced statistics, instead its mean 
value is used to scale the extracted asymmetries A. In order to avoid false asymmetries, care was taken to 
select only such data for the analysis for which the spectrometer performance was stable in consecutive 
periods of data taking. This was ensured by extensive data quality tests. In a first step, the detector 
performance was investigated on the time scale of a SPS extraction, typically 4.8 s every 16.8 s. Quan- 
tities directly linked to detection and reconstruction stability were studied, like number of interaction 
vertices, number of tracks, number of clusters in the calorimeters and trigger rates. Time intervals with 
irregularities in these variables were discarded from further analysis. In a second step, the distributions 
of DIS and SIDIS variables were investigated. In order to ensure sufficient statistics, this was done on 
the basis of a run that typically consists of 200 beam extractions. Runs showing distributions that are 
statistically incompatible to the majority of runs were also discarded from the analysis. The remaining 
data sample was carefully scrutinised for a possible systematic bias in the final asymmetries. Here, the 
two main sources for uncertainties are false asymmetries, which can be evaluated by combining data 
samples with same target spin orientation, and effects of the acceptance, which can be evaluated by com- 
paring subsamples corresponding to different ranges in the azimuthal angle of the scattered muon. The 
resulting systematic uncertainty for each data point amounts to about 75 % of the statistical error for both 
targets. An additional scale uncertainty of 5.4 % for the ^LiD and 2.2 % for the NH3 target accounts for 
uncertainties in the determination of target polarisation and target dilution factor ll36l . 



9 




X z M,, [GeV/c2] 

Fig. 4: Deuteron and proton asymmetries, integrated over the angle 0, as a function of x, z. and M/,/,, for 
the data taken with the ^LiD (top) and NH3 target (bottom), respectively. The open data points in both 
asymmetry distributions vs. M/,/, include all hadron pairs with an invariant mass of Mhh > 1-5 GeV/c^. 
These pairs are discarded for the two other distributions, which are integrated over M/,/,. The grey bands 
indicate the systematic uncertainties, where the last bin in Mhh is not fully shown. The curves show the 
comparison of the extracted asymmetries to predictions ll37l |38]| made using the transversity functions 
extracted in Ref. ifTSl (solid lines) or a pQCD based counting rule analysis (dotted lines)Pl 



5 Discussion of Results 

The resulting asymmetries are shown in Fig. [i] as a function of x, z and M/,/, for the ^LiD (top) and NH3 
(bottom) targets, respectively. For ^LiD, no significant asymmetry is observed in any variable. For NH3, 
large negative asymmetries are observed in the region x > 0.03, which implies that both transversity 
distributions and polarised two-hadron interference fragmentation functions do not vanish. For x < 0.03, 
the asymmetries are compatible with zero. Over the measured range of the invariant mass Mhh and z, the 
asymmetry is negative and shows no strong dependence on these variables. 

When comparing the results on the NH3 target to the published HERMES results on a transversely po- 
larised proton target ll28l . the larger kinematic region in x and Mhh is evident. However, both results can- 
not be directly compared for several reasons: (1) The opposite sign is due to the fact that in the extraction 
of the asymmetries the phase 7i in the angle <pRs is used in the COMPASS analysis; (2) COMR\SS calcu- 
lates asymmetries in the photon-nucleon system, while HERMES published them in the lepton-nucleon 
system; both agree reasonably well when including D„„ corrections for HERMES; (3) HERMES uses 
identified tt+tt^ pairs and COMR^SS h^h^ pairs; (4) COMPASS applies a minimum cut on z, removing 
a possible dilution due to contributions from target fragmentation. 

A naive interpretation of our data, based on Eq. Q and on isospin symmetry and charge conjugation, 
yields Di „ = Di ^ and H^^^ = —Hf^ 1271 . When considering only valence quarks, the asymmetry A^^jf 

is proportional to [/z -t- h^^H^^ for the deuteron target, while for the proton target Ajj^^^^ cx; 

Therefore, like in the case of the Collins asymmetry, the small asymmetries observed for the deuteron 
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target imply h" « — /zp while the sizeable asymmetries for the proton tai^get imply a non-vanishing u 
quark trans versity h\. 

In an early theoretical approach, a strong s and p wave interference in two-pion production, known from 
two-meson phase shift analyses, was considered for a transversely polarised proton target [21 J. As a 
result, a sign change of the two-pion asymmetries in the vicinity of the p mass was predicted. From our 
data, however, such a sign change can clearly be excluded. 

In a different approach ||25] - |27]| . all two-hadron fragmentation functions for two-pion production were 
calculated in the framework of a spectator model for the fragmentation process. In the mass range up 
to M/if, = 1.3 GeV/c^, all relevant two-pion channels with relative s and p waves were considered and 
the parameters of the model were tuned to fit the output of the PYTHIA [39] event generator tuned 
to HERMES kinematics. Predictions were made for the IFF as well as for Di „/ and Diji and in 
Ref. ilJi the expected asymmetries for COMPASS using deuteron and proton targets were calculated as- 
suming different models for the transversity distributions. Very recently [37], these parametrisations of 
the two-hadron fragmentation functions were used together with the transversity distributions extracted 
from single-hadron production ifTSll to make predictions for both proton and deuteron targets using the 
kinematic range covered by COMPASS. The applied cuts of 0.2 < z < 0.9, 0.3 < M/,/, < 1.6 GeV/c^, 
0.003 < X < 0.4, and 0.1 < y < 0.9 assure a wide overlap with the presented data. The solid lines in 
Fig. [4] show the good agreement of these predictions in comparison to the data obtained using the ^LiD 
and NH3 targets. For ^LiD, the predicted asymmetries vs. z and M/,/, are smaller than 0.01, in good 
agreement with the data. Only for x > O.l the predicted asymmetries start to exceed 0.01, a trend which 
can neither be confirmed nor excluded by the data. For NH3, significant asymmetries are predicted and 
especially the shape depending on x is well described. Also the dependences on z and M/,/, are reasonably 
described, showing a rather weak dependence on z as well as the enhanced asymmetries close to the p 
mass. The predicted kink at the p mass is not visible in the data but can also not be excluded. 
The same parameterisations of the two-hadron FF were also used in Ref. f38l, together with two models 
for the transversity distributions, a SU(6) quark-diquark model and a pQCD-based counting rule anal- 
ysis. For the latter, the results using proton or deuteron targets are shown by the dotted lines in Fig. [4] 
The cuts applied in Ref. |[38l were similar to those used in Ref. ||37]| . except that the invariant mass range 
was restricted to M/,/, < 1 .0 GeV/c^, the z region was extended to 0. 1 < z < 0.9 and no explicit x cut was 
applied. The resulting predictions describe the data fairly well. While for the ^LiD target the predicted 
asymmetries are still in agreement with the data, for the NII3 target they tend to slightly overestimate the 
measured asymmetries. 

6 Conclusions 

In this Letter we present the results on azimuthal asymmetries measured in two-hadron production in 
semi-inclusive deep inelastic scattering using transversely polarised ^LiD (deuteron) and NH3 (proton) 
targets. These asymmetries provide independent access to the transversity distribution functions com- 
plementary to that using the Collins effect in single-hadron production. For the deuteron target, no 
significant asymmetries are observed. For the proton target, sizeable asymmetries are measured in all 
three kinematic variables, x,z and M/,/,, indicating non- vanishing M-quark transversity and two-hadron 
interference fragmentation functions. No dependence of the asymmetries is observed for the variables z 
and M/,/,, on which the interference fragmentation function directly depends. Especially, the change 
of sign in the vicinity of the p mass predicted by Ref. 11211 is excluded. The dependence of the asymmetry 
on the Bjorken variable x, observed on the proton target, constitutes independent experimental informa- 
tion on the M-quark transversity distribution function h"{x). In conjunction with the unique results from 
the deuteron target, also the J-quark transversity function hy{x) can be accessed. In an interpretation 
based on valence quarks only, the proton data indicate the M-quark transversity function to be sizeable, 

^The sign of the original predictions was changed to accomodate the phase K in the definition of the angle used in the 
COMPASS analysis. 
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while the deuteron data indicate an approximate cancellation of the u- and J-quark transversity functions. 
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